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ABSTRACT 

Based on spectroscopic signatures, about one-third of known H2O maser sources in active 
galactic nuclei (AGN) are believed to arise in highly inclined accretion disks around central 
engines. These "disk maser candidates" are of interest primarily because angular structure and 
rotation curves can be resolved with interferometers, enabling dynamical study. We identify five 
new disk maser candidates in studies with the Green Bank Telescope, bringing the total number 
published to 30. We discovered two (NGC 1320, NGC 17) in a survey of 40 inclined active 
galaxies (v sys < 20000 km s _1 ). The remaining three disk maser candidates were identified in 
monitoring of known sources: NGC 449, NGC 2979, NGC 3735. We also confirm a previously 
marginal case in UGC 4203. For the disk maser candidates reported here, inferred rotation 
speeds are 130-500 km s _1 . Monitoring of three more rapidly rotating candidate disks (CG 211, 
NGC 6264, VV 340A) has enabled measurement of likely orbital centripetal acceleration, and 
estimation of central masses (2-7xlO 7 M ) and mean disk radii (0.2-0.4 pc). Accelerations 
may ultimately permit estimation of distances when combined with interferometer data. This 
is notable because the three AGN are relatively distant (10000 < v sys < 15000 km s" 1 ), and 
fractional error in a derived Hubble constant, due to peculiar motion of the galaxies, would 
be small. As signposts of highly inclined geometries at galactocentric radii of ~ 0.1-1 pc, disk 
masers also provide robust orientation references that allow analysis of (mis) alignment between 
AGN and surrounding galactic stellar disks, even without extensive interferometric mapping. 
We find no preference among published disk maser candidates to lie in high- inclination galaxies. 
This provides independent support for conclusions that in late-type galaxies, central engine 
accretion disks and galactic plane orientations are not correlated. 

Subject headings: galaxies: active — galaxies: Seyfert — ISM: molecules — masers 



1. INTRODUCTION 

Extragalactic H2O maser emission from active galactic nuclei (AGN) in spiral galaxies in some cases 
traces highly inclined disk structures at radii ~ 0.1-1 pc from the central engines. Cases est ablished by di- 



rect interferometric imaging that resolves position- velocity structure include NGC 4258 (e.g.. lMiyoshi et al 



1 Current address: Goldman Sachs, London, UK. 
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19951 : lArgon et alJl2007h . NGC 1068 dGreenhill & Gwinnlll997i). NGC 3079 (Yama uchi et al. 2004; Kon- 



dratko, Gr eenhill, fc Moran 20 05), Circinus (IGreenhill et al. M2003eJ >. UGC3789 (jReid et al J 12009 ). and 
NGC 6323 (jBraatz et al. 1 120071 ). These "disk masers" have common spectroscopic characteristics, highly 
red- and blueshifted emission lines, called "high- velocity" emission lines, that more or less evenly bracket 
the systemic velocity, and emission close to the systemic velocity, called "low- velocity" emission lines. Both 
are defined in the frame of the rotating disk, projected along the line of sight. We use this terminology 
of "high-" and "low-" velocity lines throughout this paper. This template is also observed in other H2O 
masers for which imaging data are not available, and it can be a pow erful diagnostic. In total, at least 30 



Greenhill. Tilak. &: Madeiski 



2008, 



published masers (out of 107 known) share these characteristics (e.g., 
and references therein). The high- velocity red- and blueshifted emission has been observed to be offset 
on the order of 10 2 -10 3 km s _1 , among sources. We note that very long baseline interferometric (VLBI) 
imaging is essential to confirm the identity of disk masers. Hence, we have used the terminology "disk 
maser candidates" for sources with spectroscopic identification only. We also note that VLBI imaging of 
some disk masers shows signs of additional maser emission arising in outflows (e.g., Circinus) and jets (e.g., 
NGC 1068). 

H2O disk masers highlight structures oriented edge-on, or nearly so, and in close proximity to central 
engines. This is plausible b ecause large r eservoirs of molecular gas along lines of sight are required to 
generate detectable emission (jElitzur 111982 ). Among all H2O masers in AGN, 95% exhibit X-ray absorption 
column densities > 10 23 cm -2 (Greenhill et al. 2008). Among disk masers, the fraction is nearly as high, 
and the distribution of obscuring column densities peaks above 10 24 cm -2 . For thin disks, maser emission is 
beamed narrowly about t he midplane, into a s olid angle <C 47r, as dictated by requisite velocity coherence 
along lines of sight (e.g., iMivoshi et al.lll995l ). Geometry limits observable emission to disks that are 
highly inclined. Although c ommon association has been well recognized with optical type-2 and partly 
obscured type-1 objects (e.g-. lBraatz. Wilson, fc Henkel Ill997bl ). disk maser emission is a sharper indicator 
of (inclination) orientation for nuclear structure than optical classifications are. In addition, it is a more 



direct indicator than, e.g., models of Fe Ko emission from hot accretion disks (INandra et al. 1119971 ). the 



differential brightness of radio jets, and models of kiloparsec-scale ionization structures. Comparisons of 
i nclination complement analyses of position angles and naturally enable analy ses of relative 3-D orientation 
(|C!arke. Kinnev. fc Pringlelll998l : IPringle et aUll99fll : iNagar fc Wilsonlll999l ). 

We conducted a survey of AGN to detect new disk masers, selecting galaxies with inclined stellar disks 
(inclination angle i > 70°). We detected two disk maser candidates and one unclassified maser (Section 3.1). 
We also conducted deep integrations to identify indicators of disk maser emission among previously known 
sources, detecting another three (Section 3.2), and monitored another three known disk maser candidates 
to measure centripetal acceleration due to disk rotation (Section 3.3). We discuss estimation of black hole 
masses and disk radii in Section 4 as well as the distribution of inclinations among disk maser candidate 
hosts and misalignment between central engines and galactic stellar disks. 
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2. OBSERVATIONS 

The survey was conducted during 2004-2005 with the Green Bank Telescope (GBT) of the NRAO0 
using two 200-MHz IFs configured to overlap by 50 MHz, in each polarization. We tuned the center of 
the resulting 350-MHz instantaneous bandwidth to the systemic velocity of each target source (350 MHz 
~ 4880 km s" 1 for a representative recessional velocity of 5000 km s^ 1 , assuming the optical definition of 
Doppler shift). The spectrometer delivered a channel spacing of 24.4kHz. The dedicated observation of 
known H2O maser sources was conducted between 2003 and 2006 with the same setup, except for galaxy 
CG211 (2003 December 08), where the two IFs were configured to overlap by 20 MHz and the IF at the 
lower velocity was centered on the systemic velocity of the galaxy. 

To obtain total-power spectra, we nodded the telescope by 3' every 2 minutes between two positions 
on the sky so that each target was always present in the on-axis or off-axis beam. System temperatures 
(T sys ) were measured against a calibrated noise source injected at the receiver and ranged from 23 to 75 K 
depending on elevation and weather conditions. To calibrate the spectra to flux density units, we used the 
gain curve obtained by the GBT staff based on measurements of opacity-corrected antenna temperature 
for NGC 7027 at ~ 1.4cm (R. Maddalena, private communication). By comparing maser line amplitudes 
among beams and polarizations, we estimate that the calibration of the T sys is accurate to within 30%, an 
uncertainty that dominates the error budget for the flux density scale (and is probably due to calibration 
uncertainties for the noise sources). In comparison, antenna pointing corrections were obtained every ~ 30 
minutes and were typically < 6", which corresponds to a loss less than 8% in source flux density for a 
36" beamwidth (FWHM) at 1.3 cm. Additional losses due to pointing errors induced by varying wind 



conditions were less than 20%, which we estimate using the error model of lCondonl (J2003J) and the bottom 
90% of measured wind speeds at GBT weather station 2, for each observation. (Our estimated losses are 
conservative, considering that wind speeds on the feed arm itself were typically lower.) 

Data were reduced using custom scripts written in the Interactive Data Language. We subtracted 
a running boxcar average of width 6.25 MHz to remove systematic baseline structure from the total- 
power spectra, excluding velocity intervals with signal that was detectable above the structure seen in raw 
integrations. The resulting la noise levels attained in an integration of ~ 30 minutes and corrected for 
atmospheric opacity estimated from tipping scans (from 0.03 to 0.2) were 4-14 mJy in a 24.4-kHz channel. 
The spectra presented here have been Hanning smoothed (multiple times) to a resolution of 108 kHz, 
though channel spacing of 24.4 kHz is shown in the figures. 



The National Radio Astronomy Observatory is operated by Associated Universities Inc., under cooperative agreement 
with the National Science Foundation. 



4 



3. RESULTS 



3.1. Survey — Three New Masers 

In a survey of 40 inclined (i > 70°) AGN with v sys < 18000 km s^ 1 selected from the NED (Ta- 
bled]), we have detected three new H2O maser sources in Seyfert 2 systems: NGC 1320, NGC 170, and 
IRAS 16288+3929 (Figure [T]). The three discoveries were subsequently confirmed with the Very Large 
Array (VLA) of the NRAO using a 6.25-MHz observing bandwidth and 97.7-kHz channels (NGC 1320, 
IRAS 16288+3929) or with the GBT on a different day but using the same spectrometer configuration 
(NGC 17, IRAS 16288+3929). Each of the masers is probably associated with nuclear activity. The 
NGC 1320 and IRAS 16288+3929 maser positions measured with the VLA lie <C lcr from the optical posi- 
tions of the nuclei (Table[T|), or less than 300 pc. The NGC 17 maser is three t o four times more distant than 



the fa rthest H2O maser known to be associated with star formation (Arp 244; iDarling. Brogan. fe Johnson 
20081 ). yet it is comparably bright. 



Spectroscopic classification of sources as disk maser candidates is most suggestive when three line 
complexes are detected, one close to the systemic velocity and two bracketing it, with similar but opposite 
velocity offsets. Approximate symmetry in this respect is indicative of a disk, and a minimum offset on 
the order o f 200 km s" 1 is reasonable to exclude misidentification of emission driven by nonnuclear star 
formation (jGreenhill 1 120071 ). However, detailed geometry of the underlying accretion disks, orientation 
with respect to the line of sight, and sensitivity limitations may result in fewer than three detectable line 
complexes. In these cases, interpretation of the spectra depends critically on estimates of systemic velocity. 

The spectrum of NGC 1320 is relatively symmetric about a mean of ~ 2790 km s" 1 (Figured]), with 
emission extending outward ~ 170 km s _1 from peaks at 2614 and 2969 k m s -1 . However, t he mean 
velocity is offset ~ 70 km s _1 from the most precise op tical systemic ve l ocity ( Huchra et al.lll993l ). though 
agreement appears to be good with the estimate of iBottinelli et al.l (|1992l ). The approximate mirror 
symmetry and breadth of the line complexes is suggestive of 150-180 km s -1 rotation velocity, though it 
would be twice as great if one of the detected line complexes were to correspond to systemic emission. Time 
monitoring might distinguish between scenarios if a third, blueshifted, variable line complex is detected or 
if secular drift in velocity is detected in one of the known complexes. 

Classification of NGC 17 is somewhat less certain because of a spread in systemic velocity estimates 
that spans the separation of the two detected line complexes ~ 300 km s" 1 (Figure [I]). However, we 
note that the most precise estimate (IRothberg <fe Joseph! 120061 ) . obtained from dynamical modeling of 
the galaxy, corresponds well to the peak of the stronger of the two line complexes. We infer from this 
circumstance that NGC 17 is likely to be a disk maser, with rotation velocity ~ 300 km s" 1 . As for 
NGC 1320, monitoring may increase certainty s hould a blueshifted line complex or secular drift be detected . 



Nonetheless, bec a use NGC 17 is a merger (e .g., IVorontsov-Vel'Yaminovlll977l ; iMazzarella &: Borosonlll993l ; 



Hunt et al.lll999l ; ISchweizer &; Seitzerl 120071 ) . interpretation of the spectrum requires some caution. The 
velocity offset b etween the two spectral f eatures is comparable to the stellar velocity dispersion, a Q = 
201 ± 8 km s" 1 (jRothberg fc Joseph) IgOOJ ), and the maser lines may arise in different compact structures 



2 DupIicat e designation, NGC 34, originates from a suspected observer error recorded in iDrever I (188l|) an d corrected by 
Howel (|l9O0h . Dual designation does n ot stem from NGC 1 7 being a merger with a complicated brightness distribution. 
However, see discussion of NGC 17/34 in Hopkins et al. ( 20081 ). 
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within the system. 



For IRAS 16288+3929, classification is problematic because only a single spectral feature is evident 
(~ 8920 km s _1 ), and it is displaced ~ 200 km s _1 from available optical systemic velocities (Figure [1]). 
The emission may represent high- velocity blueshifted emission or activity unrelated to any accretion disk. 
An alternative origin could be jet activity driven by the central engine. We note comparable v e locity 



offsets for known "j e t masers" (NGC 1068, NGC 1052, Mrk348, and M 51 ; lOallimore et alJll99fl . 



Claussen et al 



1998 



Kameno et al. 


2005: 


Peck et al. 


2003: 


Hagiwara 


2007) 



2001 



cases, spectral features are characteristically at least 50 km s 1 wide, at least an order of magnitude broader 
than the linewidth for the known maser feature in IRAS 16288+3929. 

NGC 1320 and NGC 17 were targeted but not detected in previous surveys for maser emission. Braatz, 
Wilson, & Henkel (1996) and Kondratko et al. (2006a) report observation of NGC 17 with la noise levels of 
73 and 16mJy in 0.66 and 1.3 km s _1 spectral channels, respectively. These correspond to signal-to-noise 
ratios (SNR) of < 1 for the observed line flux of llmJy (peak) over 1.5 km s _1 . In the case of NGC 1320, 
these earlier studies obtained la noise levels of 61 and 19mJy in 0.84 and 1.3 kms -1 spectral channels, 
respectively. The SNR for the observed line flux of 37mJy, at 1.5 km s" 1 resolution, would have been less 
than two. 



3.2. Three New Disk Maser Candidates (NGC 449, NGC 2979, and NGC 3735) among Known Sources 



Through deep integration, we detected new features in the spectra of masers not previously recognized 
to have high- velocity emission and evidence for origins in accretion disks (Figures[2l O Table[2|) . Noise 
levels were 2-4 mJy (la) referenced to outside the atmosphere. We comment on individual systems below. 
The first three we discuss are the new disk maser candidates. 



NGC 449 (Mrkl) (Figure [2]). We detected an emission complex at ~ 4860 kms" 



a ~ 



20mJy com- 



ponent at ~ 4600 kms" 1 , and a new weak line complex (7mJy peak in a 1.5 kms" 1 spectral channel) at 
~ 4720 k m s" 1 . The complex at ~ 4860 km s -1 h as been present in the spectrum of the source since its 
discovery ( Braatz et al. 1996 : iBraatz et al.ll2003l ) while the component at ~ 4600 kms" 1 was marginally 



detected first by 



Braatz et al. 



(|2003l ). The newly detected emission is approximately halfway in velocity 
between the two previously known spectral complexes. Although the position of the maser emission with 
respect to the nucleus has not been measured, the spectrum shows the archetypical spectral signature of 
emission from a highly inclined accretion disk. In particular, considering the uncertainty on the systemic 
velocity of the galaxy (Figure[2j Table[2]), we associate the complex near 4720 km s _1 with low- velocity 
disk emission and infer an orbital velocity (which we take to be the maximum velocity with respect to 
the systemic velocity) of ~ 150 km s _1 . It is notable that the suggested systemic velocity of the galaxy 
would be blueshifted by a few times 10 km s" 1 from optical measures and ~ 100 km s _1 from the estimate 
obtained from HI data. 

NGC 2979 (Figure [2]). Emission was newly detected at a level of ~ 3mJy at 2390 km s _1 (cf. 
Greenhill et al.1 l2003bl ). which is blueshifted by ~ 300 km s" 1 from estimates of the systemic velocity. 
The new spectrum is approximately centered on the optical systemic velocity estimate of 2720 + 15 km s" 1 
( Fisher et al,lll995l ). which is suggestive of emission from a highly inclined disk with ~ 300 km s" 1 orbital 
velocity, though extended line wings may suggest a somewhat more rapid rotation. This interpretation is 
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consist ent with the known coincidence of the maser source position and the active nucleus (jGreenhill et al 
2003bl ), Asymmetry in the maser spectrum is pronounced but difficult to explain without imaging data 
that would enable geometric modeling of the disk. However, we note that the archetypal disk maser in 
NGC4258 exhibits about a 10 times imbalance in the flux densities of red- and blueshifted emission. In 
the new spectrum, the velocity range of emission > 10 mJy is consisten t with a weak plateau observed in 
the 6 times lower sensitivity discovery spectrum (jGreenhill et al.ll2003bl ). 

NGC3735 (Figure EJ). W e detected narrow lines at 2444, 2484 2784 , and 2896 km s" 1 that were not 
observed in previous spectra (jGreenhill et al.l Il997bl : iBraatz et al.1 120031 ) . An additional complex of weak 
lines may be present between 2 830 and 2870 km s -1 . Doppler shifts from systemic velocities derived from 
HI data are ± 200-250 km s _1 (jStavelev-Smith &: Davies 1988). This bracketi ng of systemic emission and 
known angular coincidence of the maser and nucleus (jGreenhill et al.lll997bl ) are suggestive of emission 
from a close-to edge-on accretion disk. 

UGC 4203 (Mrk 1210) (FigureEJ)- We detected line complexes bracketing the systemic velocity, at 3755- 
3861 km s" 1 and 4190-4350 km s _1 (Figure[3j), as well as isolated narrow lines distributed between these, at 
4106 km s~ 1 (~ 70mJy), 4134 km s _1 and 3992 kms -1 (5-10 mJy ). The peaks of the red- and blueshifte d 



complexes were detected with lower signal-to-noise ratio previously (jBraatz et al. 1119961 : IBraatz et al,ll2003l ) , 
though the symmetry of the spectrum around the systemic velocity (Table[2]) and suggestion of emission 
from an inclined disk were less pronounced. The inferred orbital velocity is ~ 300 km s -1 , depending on 
the adopted systemic velocity. 

Arp299 (NGC3690) (Figure EJ). The spectrum contains a broad feature with FWHM ~ 230 kms" 1 
and overlying narrow Doppler comp onents that were n ot evident in the previous spectra of the source 
obtained with 4.3 km s _1 channels ( Henkel et al.l 120051 ) . Maser emission is known to arise from three 
regions in t his merging system , the nuclei of NGC 3690 and IC 694 and an internuclear radio continuum 
"hotspot" (jTarchi et al.l 120071 ). Both nuclear maser sources are ~ 200 km s _1 wide, which is notable 



because large widths have been associated at least circumstantially with jet activity (see lTarchi et al.ll2007l ). 
as exemplified in three AGN where maser and radio continuu m emission have been mapped relative to one 



another inter ferometrically: NGC 105 2, NGC 1068, Mrk 348 (jClaussen et al 



Peck et al 



2003 



1998 



Gallimore et al 



2001 



Kameno et al.ll2005l . also review by Greenhill 2007). Notable time variability has been 



observed toward these masers, and position-resolved monitoring of Arp299 may offer useful diagnostics. 
NGC 4293 (Figure[3|). The spectrum exhibits emission distributed over ~ 29 km s _1 but that unevenly 



brackets the most precise estimate of systemic velocity (jdi Nella et al.1 Il995l ) . An additional emission 
feature at ~ 1200 km s~ 1 may be present at the ~ 3a level. The case for origin in a disk is marginal due 
to uncertainty in the full velocity range of emission and the degree of symmetry. 

NGC 0235 A (Figure H]). Dominated by a single spectral feature redshifted by ~ 300 km s" 1 from the 
systemic velocity, interpretation of the spectrum is difficult. Goals of further monitoring should include a 
search for blueshifted emission. 

NGC 2824 (Figured]). As observed in the case of Arp 299, the NGC 2824 maser exhibits narrow Doppler 
componen ts overlying a broad c omponent (FWHM ~ 150 kms -1 ). This is consistent with the discovery 
spectrum (jGreenhill et al.ll2003bl ). though the newly recognized zero-intensity widt h is substantially broad - 
ened, ~ 450 km s _1 . The maser emission is associated in position with the AGN ( Greenhill et al.ll2003bl ). 
and the spectrum is suggestive of jet activity, but the interpretation is uncertain. 
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IRASF 19370-0131 (Figur e HI) . Observed emiss ion is limited to a range ~ 50 km s _1 centered on the 
systemic velocity estimated bv lStrauss et al.l (|1992l ). Detection of secular drift in feature velocities, as well 
as emission over a broader range would be useful goals for future monitoring. 



3.3. New Centripetal Accelerations for Three Disks 

Deep integrations were also used in monitoring three masers reported by Kondratko, Greenhill, & 
Moran (2006b), which display among the largest Doppler shifts for high-velocity emission: CG 211, NGC 6264, 
and VV340A (Tabled Figure[5|). We detected secular drifts in low- velocity features, which probably reflect 
centripetal acceleration due to disk ro tation (Table [H Figs. [6] a nd [T]). For each source, we used an iterative 



least-squares technique, adopted from iHumphrevs et al.l (|2008l ). to decompose the time-series spectra into 



Gaussian components, solving simultaneously for the light curve, time-varying line width, velocity at a 
reference epoch, and velocity drift of each component. (Each component was assumed to have its own 
fixed drift rate.) We constrained the widths of most Gaussian components to be < 5 km s" 1 , motivated 
by empirical observation. A first-order polynomial and a wide (~ 25kms~ 1 ) stationary Gaussian com- 
ponent were used to model the broad low- velocity plateaus in VV340A and CG211, respectively. The 
fitted accelerations were not sensitive to inclusion of the broad components, though reduced x 2 estimates 
were improved. We obtained drifts of 0.1-1.9 km s~ 1 yr~ 1 with uncertainties of 0.005-0.09 km s yr _1 , 
1.2-2.1 km s _1 yr _1 with uncertainties of 0.1-0.6 km s _1 yr _1 , and 3.8 =L 0.1 km s _1 yr _1 for CG211, 
NGC 6264, and VV340A, respectively (Table 0}. 



4. DISCUSSION 

4.1. Black Hole Masses, Disk Radii, and Their Implications 

Mean measured centripetal acceleration, a, combined with mean orbital velocity inferred from single- 
dish spectra, v, can be used to estimate central-engine gravitational masses and accretion disk radii. The 
mean acceleration among low- velocity features may be expressed as a/ = vf/Ri, where v\ is an effective 
mean rotational velocity and R[ is the effective mean radius, both for material contributing low-velocity 
emission. The mean rotational velocity inferred from spectra may be expressed as v\ = GM/ Rh, where G is 
the gravitational constant, M is black hole mass, and Rh is the effect mean radius for material contributing 
high- velocity emission. For edge-on orientation, we obtain M = vf l vf(Rh/Ri)/Gai. Unfortunately, without 
VLBI, we cannot determine Rh and Ri independently or, therefore, Vh and vi independently. Hence, for the 
purpose of calculation, we assume that the low- and high-velocity emission arise at the same effective mean 
radii {Rh ~ Ri = R) and that the effective mean rotation velocity of accelerating low-velocity material 
is the same as the measured mean rotational velocity, Vh ~ V\ = v. We then obtain M = v 4 /Gai and 
R = v 2 1 a. These values are reported in Table [5j VV340A may have the largest central mass among known 
maser galaxies. 

We note that the scatter in the velocity drift measurements for CG 211 is substantially larger than the 
individual formal uncertainties. Since a = GM/r 2 for low- velocity material and hence Aa/a = — 2Ar/r, 
the scatter in the measured accelerations of Aa/a ~ 1.0 might be due to low-velocity emission that 
populates a wide range of radii within the accretion disk, i.e., Ar/r ~ 0.5. In our simple model, the radial 



s 



spread among high- velocity material might be expected to be Av/v = 0.5 Ar/r ~ 0.25. The observed value 
of Av/v in the high- velocity components is about 0.3 (see Figure 5). This suggests that the radial spread 
is about the same for both high- and low-velocity masers. For NGC6264, Aa/a ~ 0.4, so Ar/r ~ 0.2, 
and Av/v for the high-velocity masers is expected to be about 0.1. However, inspection of the spectrum 
(Figure [5]) suggests that Av/v ~ 0.5. Hence, the high- velocity masers probably occupy a larger range of 
radii than the low-velocity masers. Note that the case of NGC 6264 may be similar to that of NGC 4258, 
where the radial range of the high-velocity features is also greater than that of the low- velocity features 
(jHerrnstein et al.ll2005l ) 

In each case where we report centripetal acceleration, it has been positive, i.e., spectral components 
drift towards larger velocities. This is consistent with the emission being from material on the near side 
of the disks. Material behind a dynamical center, which would exhibit negative centripetal accelerations, 
has not been detected toward any known disk maser or disk maser candidate. Although the relatively 
low signal-to-noise achieved for CG211, NGC 6264, and VV340A would make recognition of "back-side" 
emission relatively difficult, we can exclude the possibility that flux densities are equivalent to the front- 
side emission. It is difficult to know whether the apparent absence of back-side emission may be due to a 
maser radiation propagation effect (e.g., competitive gain diminishes inward-directed beams) or differential 
absorption (e.g., inward-directed beams pass through regions of enhanced scattering or absorption in the 
vicinity of a central engine; inward-directed beams cross the full breadth of a disk and are attenuated by 
high emission-measure material at large disk scale heights, as in Herrnstein, Greenhill, & Moran 1996). 

The three masers in Table[5] are priority targets for interferometric programs that are intended to 
estimate geometric distances (e.g., Herrnstein et al. 1999) and (in combination with distances for other 
maser targets) a combined Hubble constant, Ho. For a sample of N anchor galaxies, uncertainty in Ho 
scales as ~ _/V~ 0,5 (for uniform uncertainties) because errors for individual distance measurements are 
largely uncorrelated; they depend principally on the distribution of Doppler components across the disk 
face and details of disk geometries, which differ from galaxy to galaxy. 

To provide substantive independent constraint on additional cosmological parameters, such as a time- 
variable equation of state (EOS) for dark energy and curvature, estimates of Ho need to achieve on the 
order of 1% uncertainty, matching the growing accuracy of other measurem e nts, e.g., Cosmic Microwave 



Background (CMB) fluctuat i ons a nd Baryon Acoustic Oscillations (|Hu M2005I : llchikawa &: Takahashi 



Oiling 



2007 



2008 



Greenhill et al.1 l2009l . and references therein). In combination with CMB data alone, a 1% 



estimate of Hp wo uld enable a 5% constraint on a presumed-constant EOS, for a flat universe (Figure 14, 
Riess et al. Il2009l ). A sample on the order of 100 disk masers would be required to achieve 1% accuracy 
in Ho, if each distance is accurate to 10%. At present, this is at least an order of magnitude beyond the 
number of disk masers apparently suitable for distance measurement. However, a three- year survey of 3000 
galaxies at the GBT, achieving RMS noise ~ 3 mJy over 1 km s^ 1 (4 sources per hour for ~ 750^), would 
res ult in the detecti o n of 1 50 masers, assuming a conservative detection rate of one-half that reported 
bv iKondratko et al.1 (12006bh for AGN at 10000-15000 km s _1 . For a disk maser fraction of one-third, 
and a further one-quarter that appear suitable for distance measurement (as suggested in the currently 
known population), it would be possible to assemble a sample on the order of 25 objects. For 10% distance 
uncertainty per maser, this would enable estimation of Ho with 2% accuracy in the near term. Prospects 
are improved by better individual distances as may be achievable with the VSOP-2 space- VLBI mission 



and larger target samples made accessible by the Square Kilometer Array (jGreenhill ll2004l ; lMorganti et al 
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200J). 



However, peculiar radial motions of galaxies increase scatter in the Hubble relation and introduce 
partial correlation among recessional velocities. The estimate of the Hubble constant from any particular 
galaxy measurement is D/v. The sample size needed to achieve a given total uncertainty grows as the 
quadrature sum of AD/D and fractional peculiar motion, v P ec/v S y S , so that the uncertainty in a sample 



SJ 2 + (AD/D) 2 } 1/2 JV-V2. if AD/D 



of N galaxies is AHq/Hq ~ [(v pec /v, 

uncertainty in peculiar motion raises the error by less than 10%, then v pec /v sys 



0.1, and we require that the 
< 0.05. Outside of clusters, 



v pec < 500 km s 1 is typical (jSpringob et al.l 120071 ). We take this to be an RMS value and infer that 
v sys > 10 4 km s _1 is needed. For m aser hosts inside ric h clusters, barycenter velocities can be estimated 



with uncertainty < 200 km s 1 (e.g.. lCrook et al. 1120071 ) and used in place of individual galaxy recessional 



velocities, though cluster peculiar motions would still deserve attention. 

It would be desirable to target relatively nearby galaxies that are in "quiet" portions of the Hubble 
flow, to reduce statistical and nonstatistical errors (e.g., UGC 3789; lMastersll2005l ; iReid et al.ll2009i ). How- 
ever, a broad distribution of sample members in velocity and angle on the sky is needed. Subtraction 
of cosmic flow models to at least partially correct for Vp ec reduces the lower limit on acceptable v S y S and 
admits galaxies in not so quiet regions. Using Tully- Fisher dis t ances to estimate peculiar velocities and 
construct a multi-attrac t or mo del inside 6000 km s , iMastersI ( 20051 ) obtained a 160 ± 20 km s" 1 RMS 
residual. lErdogdu et al.1 (|2006l ) report a somewhat less direct error estimate, a scatter of 100-150 km s _1 
around a flow model derived from galaxy recessional velocities (as opposed to peculiar velocities), up to 
~ 16000 km s . Adopting a characteristic 160 km s _1 uncertainty, the lower limit in v sys is 3200 km s , 
for maser distances good to 10% and the nominal limit of ~ 10% contribution by peculiar motions to the 
error budget. The technique works for maser hosts away from the Great Attractor and the galactic plane 
(|6| > 20°; Masters 2005), where flow models are not well constrained. 



4.2. Accretion Disk and Host Galaxy Misalignments 



In studies of the relative orientation of AGN central engines and host galaxies, radio jet position angles 
at parsec and kiloparsec scales are used to infer the orientations of the accretion disks that launch them. 
Ionization cone position angles are used to infer the orientations of larger scale, dusty, geometrically thick 
structures at larger r adii. Comparisons of posit i on angles suggest that accretion disks are not aligned with 



galactic spiral disks (lUlvestad fc Wilson! Il984j : iNagar fc Wilson 



1999 



Kinney et al. 



2000 



Schmitt et al 



20031 : iMiddelberg et al.l 12004 ) . At the same ti me, the accretion dis ks appear reasonably well aligned with 
the structures that collimate ionization cones (jSchmitt et al.ll2003l ). though the sample is small. 

Because interferometric mapping has demonstrated that disk maser emission marks molecular struc- 
tures that are oriented close to edge-on, the emission provides a relatively precise indicator of orientation, 
via inclination rather than position angle, as is otherwise the case principally for radio jets and ionization 
cones. Thus, estimation of stellar disk inclination in spiral galaxies that host disk masers enables im- 
mediate estimation of (mis) alignment between parsec-scale and kiloparsec-scale structures. The approach 
complements other techniques, in terms of orientation angle (position angle vs inclination) and of radius, 
where the range of radii for maser emission lies outside the zone where accretion disks act to collimate radio 
jets. Perhaps the fullest case studies would be possible for galaxies with observable radio jets, disk maser 
emission, well-resolved ionization cones, accurate B-band isophotal radii (for detection of spiral arms), 
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and a reasonably firm morphological type. Among nearby galaxies, inclinations could be obtained from 
detai led modeling of spectroscopic data for well-resolved stellar disks and circumnuclear molecular rings 
(e.g JCurranlhoOoh . 

A suitable sample for study may be compiled directly from lists of known disk maser candidate hosts 
that are spirals. The spectroscopic identification of masers that trace underlying disk structure is sufficient. 
Galactic inclinations, to first order, may be inferred directly from axial ratios, i = cos -1 (6/a), though 
detailed dynamical modeling would be preferred to obtain the most accurate values. 

Assuming the spectroscopic identifications are accurate, then there is no apparent preponderance of 
edge-on galactic disks among disk maser candidate hosts. For inclinations < 70-80°, the distribution of ori- 
entations for the disk maser candidate matches what would be anticipated for random galaxy orientations 
(Figure [8j). To make the comparison, we assembled a sample of axial ratios for published disk maser candi- 
dates, defined with respect to surface brightness of B = 25.0 per sec -2 , for SO and la ter systems (Table [6l). 



Axia l ratios were taken from the Third Reference Catalogue of Bright Galaxies (Ide Vaucouleurs et al 



199l|). 



The apparent undercount among high-inclination systems in FigureO is probably due to overestima- 
tion of disk minor axes due to galactic bulges, whic h is anticipat e d to be greatest for nearly edge-on 



massive systems (i.e., morphologies from SO to Sa). iFouque et al.1 (|1990l ) describe an empirical correc- 
tion that is dependent on galaxy morphology, quantified in the T parameter (e.g., Ide Vaucouleurs et al 
Il99ll ). Determination of morphologies is strongly dependent on image quality and may be affected by 
subjective bias (e.g., see methodology discussion of lFukugita et al.l 120071 ). Corrected inclinations are par- 
ticularly sensitive to uncertainty in T for highly inclined systems, where forbidde n values of cosz may 
result. We estimate corrected inclinations using catalogued mor phological types (jde Vaucouleurs et al 



199ll ; iMalkan. Gor jian. fc Tarn I Il998l ; lLauberts &: Valentiin 1 1 19891 ) . For galaxies with forbidden values of 



cosi (NGC 1386, NGC4388, NGC3735), we increment T by 0.6, which is approximately the catalogued 
uncertainty for these systems. We obtain inclinations that are within a few degrees of 90°. (In light of 
the chosen coarse binning, the precise values are not critical.) As a result, agreement with what is ex- 
pected for random galactic disk orientations is obtained over the whole range (FigureEJ) . For an intensively 
studied sample that includes well-characterized inclination uncertainties for each object, formal statistical 
comparison would be possible using a Monte Carlo analysis. 

No correlation of misalignment with morphological type or mass is apparent. Disk maser candi- 
date hosts span a range of almost two orders of magnitude in black hole mass and a broad range of 
apparent morphological type, from SO (T = —2) to Sc (T < 5), excluding the merger system NGC 17. 
The median type is 2.7, with 1 < T < 4 for the middl e quartiles (Table El). Masses are inferre d 



from VLB! or acceleration measure ments: NGC 4258 (e.g., Mivoshi et al. 19951 ; Herrnstein et al. 2005 ) 



NGC 1068 (Greenhill fc Gwinnlll997), Circinus (Greenhill et al. 2003a), NGC 3079 (Yamauchi et al 



2004; 



Kondratko et alJbood). UGC3789 faeid et alE oogh. NGC 6323 feraatz et al. IhoOTl ). IC 2560 (jlshihara et al 



2001 



), Mrkl419 (|Henkel et al.1 12002 ). and CG211, NGC 6264, and VV340A (this work). 

Comparison of galactic plane misalignments inferred for disk masers and for accretion disks at smaller 
radii thus far is statistical. Direct comparison would require disk masers that have been mapped with VLBI, 
but that sample is comparatively small. In several examples, the position angles for jets (or lobes) and ion- 
ization cones are well aligned with dis k maser rotation axes, projected on the sky: NGC 3393 (Kondratko, 



Greenhill, & Moran 2008), NGC 4258 (ICecil et al.ll2000l ). and NGC 4945 (jCreenhill. Moran. fc Herrnstein 
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199 7cJ). However, misal ignment between a jet and disk maser has also been observed in the case of NGC 3079 



(jKondratko et al.ll2005l ). a nd skewing of ionization cone orientation by shadowing from a warped disk has 
been observed in Circinus (jGreenhill et al. Il2003al ). both at radii < lpc. Misalignment on the order of 30° 
is also observed in NGC 1068, between the position angle of the disk maser and th e major axis of ionized 
material just inside the inner radius for molecular emission (jGallimore et al.l 12004 ), If these examples are 
re presentative, counterp arts may contribute to scatter in the alignment of jets and ionization cones noted 
by lSchmitt et al.1 (120031 ) but ascribed to signal-to -noise and project i on eff ects. Warping as observed among 
disk masers and larger nuclear structures (e.g., ISchinnerer et al. I l2000l ) may also be responsible, where 
fueling of central engines does not depend on conveyance of material along galactic planes to ever smaller 
radii and vertical heights until it is concentrated into a coplanar accretion disk. 



5. CONCLUSIONS 

In a survey with the GBT of 40 inclined AGN, we have detected three new water maser sources 
(NGC 1320, NGC 17, and IRAS 16288-3929; Figure QJ. We classify two of these (NGC 1320 and NGC 17) 
as disk maser candidates based on their spectra. Three more cases of disk maser candidates were dis- 
covered through deep integrations toward known masers (NGC 449, NGC 2979, and NGC 3735; Figure [2]). 
One previously marginal case was confirmed (UGC4203; Figure [3]). Inferred rotation speeds are ~ 100- 
300 km s _1 . At present, the details of interpretation are sensitive to estimates of systemic velocity at 
optical wavelengths and 21 cm (the HI line). These show significant scatter for some galaxies, and in 
general, comprehensive study is mandatory for each new maser detection to resolve uncertainties in clas- 
sification. Interferometric follow-up of discoveries reported here will enable confirmation studies of disk 
geometry (e.g., radii, warping) and estimation of central engine mass, at least for the sources with large 
numbers of maser Doppler components. In general, systematic velocities that result from disk models may 
provide the best estimates of systemic velocity, and these might usefully be compared to the results of 
pointed optical and molecular line studies that include dynamical modeling. 

For three additional maser sources exhibiting high-velocity maser emission (CG 211, NGC 6264, VV340A), 
we report secular drift in multiple low-velocity features — a probable manifestation of centripetal acceler- 
ation due to disk rotation. Combined with rotation velocities inferred from spectra, these centripetal 
accelerations suggest central engine masses of a few times 10 M Q and mean disk radii of a few times 
10 _1 pc. Measurable accelerations of ~ 1-3 km s _1 yr _1 and large recessional velocities (> 10 4 km s _1 ) 
make these three sources particularly attractive for interferometric study in pursuit of robust estimates of 
geometric distances and Ho- 

As beacons of highly inclined structures in AGN, at radii of ~ 0.1-1 pc, disk masers also highlight 
those galactic nuclei in which we can define a reference plane with high accuracy and study (mis) alignment 
with surrounding kiloparsec-scale stellar disks. We find no preponderance of edge-on galactic disks among 
published disk maser candidate hosts. The comparison is in agreement with studies that infer orientations 
of accretion disks and geometrically thick tori from radio jets and ionization cones, though the maser data 
refer to a unique range of radii and probe inclination rather than position angle on the sky. 

Our studies provide more candidates for the quest to measure Ho with disk masers. The best available 
estimate of Hq is uncertain by 5% and achieved with a composite anal ysis of Cepheid a nd supernova 



data, with a zero point provided by the disk maser distance to NGC 4258 (jRiess et al. M2009I ) . Ultimately 
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accuracy on the order of 1% would match accuracies anticipated for CMB and other data sets in the near 
future and add a new independent constraint in cosmological parameter estimation. However, achieving 
this accuracy for Ho will be challenging, probably requiring detailed study of a disk maser sample at least 
three times larger than what is available today, assuming maser distances with 10% accuracy as the norm. 
The minimum sample recessional velocity would be 10 4 km s" 1 , though objects as nearby as 3200 km s" 1 
may be useful where peculiar motions are known with uncertainties as small as 160 km s~ 1 . Difficulty in 
accurately estimating peculiar motions in some volumes may restrict this extended sample, e.g., near the 
Great Attractor or where uncertainty in Tully-Fisher or Fundamental Plane distances are large. 

The graduate dissertation by PK contributed to this work. We thank M. Elvis, E. Humphreys, and R. 
Narayan for helpful discussions. This work benefitted from expert analysis of peculiar velocities and model 
flow fields offered by K. Masters. We thank M. Reid for discussions and the code used to fit the velocity 
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as overall encouragement. We thank C. Bignell for telescope scheduling that enabled the dissertation work 
and R. Maddalena for in-depth discussion of calibration. This research has made extensive use of the 
NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory (JPL). 
This work was supported by GBT student support program, grants GSSP004-0005 and GSSP004-0011, 
and in part by grant NNG05GK24G from NASA. The results presented here were obtained through NRAO 
observing programs GBT04C-031, GBT05A-015, GBT06A-056 (PI: Kondratko), and AK629 at the VLA. 
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Fig. 1— Spectra of newly discovered emission in NGC 1320, NGC 17, and IRAS 16288+3929. After 
smoothing, the spectral resolution is 108 kHz (1.5 km s _1 ), for which la noise levels are 2.8, 1.5, and 
1.3 mJy, respectively. Vertical bars indicate systemic velocity estimates. Where there is more than one, 
the thicker one represents the adopted value, selected based on c onsideration of error budget. Velocity 
estimates for NGC 1320 a re 2716 ± 29 (optical; iHuchra et aflll993h . which is adopted here, 2663 ± 16 (HI; 
Davoust fc Continil[20o3 v and 2786 ± 90 (optical; iBottinelli etal1ll992l ). (We note duplication of NGC 
designation, and the swapping of MCG identifier and velocity between N GC 1320 and NGC 1321 en tries in 
da Costa et al. 1998.) For NGC 17, the velocities are 5 881 ± 2 (optical; iRothberg fc Joseph! bood), which 
is adopted here, 5 931 ± 10 (HI: IBottinelli et al.l fl99ol) . 5772 ± 2 5 (optical; Ide Vaucouleurs et alJ fl99ll). 
5726 ± 30 (optical; bsterbrock fc DaharJll98ah. 5821 + 38 f optical: buchra et al]ll993l ), 5821 ± 44 (optical; 



da C osta et al. 1998), 5870 ± 15 (optical; ISchweizer fc Seitzerl 120071 ). and 5839 ± 80 (optical; Ijones et al 



20051 ). For IRAS 16288+3929 optical velocitie s are 9 123 ± 9 (emission lines; adopted here) and 9136 ± 27 
(absorption lin es), both from kbazaiian et al.1 (|2005l ). and 9091 ± 29 (|Zabludoff. Huchra. fc Geller lll99ol ) 
and 9161 ± 21 (|Rines et alibood ). 
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Fig. 2.— Spectra of known masers in NGC 449 (Mrk 1), NGC 2979, and NGC 3735, obtained with the GBT. 
After iterative Hanning smoothing, the effective spectral resolution is 1.5 km s -1 . Vertical bars indicate 
systemic velocity estimates; where there are more than one, the thicker one represents the adopted value 
(see Table [2]). The vertical scale was expanded to make apparent weak emission. The peak flux densities 
are noted below the vertical arrows. 
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Fig. 3.— Spectra of known masers in UGC4203, Arp299, and NGC4293, a s in Figure H Fo r Arp 299, 
systemic velocities are shown for each nucleus, NGC 3690 and IC 694. See iHibbard &; Yunl (|1999l ) for 
discussion of nomenclature and identification. 
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Fig. 4.— Spectra of known masers in NGC 0235A (ESO474-G016), NGC 2824, and IRASF 19370-0131, as 
in FigureEJ NGC0235A is part of a close pair (~ 11") with ESO474-G017, identified as a Seyfert2 object 
in the NED. The systemic velocity 6772 ± SGkms^ 1 (da Costa et al. 1998) is also close to the maser line 
velocity. 
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Fig. 5. — Variance-weighted average spectra for CG211, NGC6264, VV340A, for the epochs listed in 
TableO Spectra have been iteratively Hanning smoothed to a resolution of 1.5 km s^ 1 , resulting in 1- 
a noise levels of 0.53, 0.42, and 0.63 mJy, respectively. The vertical axes are split to emphasize Doppler 
components below 10 mJy. Vertical bars indicate systemic velocity estimates; where t here is more than one, 
the t hicker one represents the adopted value. CG 211: 14810 ± 120 km s -1 (optica l; IWeistrop &: Downes 



199ll ): NGC6264: 10141 ± 10 km s" 1 (optical; lAdelman-McCarthv et al. 120071 ): 10177 ± 28 km s 
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Falco et al. 
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( optical: beers et alJ [l99^) ■ 10161±76 km s" 1 (optical: koranvi fc Gelleribood '). 10151±45 kms" 1 (optical; 
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Fig. 6. — Spectra of low-velocity features in CG211 obtained with the Green Bank Telescope and the 
results of Gaussian component decomposition (also see Tabled]). A wide (~ 25kms~ 1 ) stationary Gaussian 
component at ~ 14951 km s _1 represents a broad low-velocity plateau present at each epoch. 
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Fig. 7. — Spectra of low-velocity features in NGC 6264 (left) and VV340A (right) obtained with the Green 
Bank Telescope and the results of Gaussian component decomposition (also see Tabled]). In the case of 
VV340A, horizontal lines to the left and right of the spectra mark the zero level, which is distinguished 
from the low-level broad emission component. 
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Fig. 8. — Histograms of inclinations for disk maser candidate host galaxies (Table |SJ), demonstrating 
misalignment with the edge-on reference planes defined by disk maser candidates on sub-parsec scales. 
(left) Inclination computed from catalogued axial ratios at B = 25 magnitude sec -2 (Ide Vaucouleurs et al 



199ll ). for published disk maser candidates that lie in spiral hosts. The dashed line represents the expec 
tation for a sample of randomly oriented galaxies, scaled to seven objects in the highest bin, (right) 



Inclin ation computed with correction for dis tortion of isophotes due to galactic bulges (IFouaue et a 
1990h. using estimates of mo rphological type T (Ide Vaucouleurs et al.lll99ll : lMalkan. Gorjian. fc Tarn I ll998l : 



Lauberts fc Valentiin Ill989l ). Shaded bars represent three galaxies with inclinations very close to 90° and 
for which a ~ lc or +0.6 adjustment to T (e.g., Sa^Sab) was required to obtain allowable values of cosi. 



Table 1. Inclined AGN Surveyed for H2O Maser Emission 



Galaxy" 


Type 


02000 
(hh mm 


b 

ss) 


#2000 6 
(dd mm ss 


) 


Vsys 

(km s" 1 ) 


Date 


rp d 

1 sys 
(K) 


\<7 e 

(mjy) 


UGC 12915 


LINER 


00 


01 


41. 


94 


+23 


29 


44. 


5 


4336 


2005- 


■02-02 


35 


5 


.6 


UGC 00050 


LINER 


00 


06 


10 


15 


+26 


09 


16 


2 


7552 


2005- 


02-01 


31 


6 


1 


NGC 17* 


Sy2 


00 


11 


06 


55 


-12 


06 


26. 


3 


5931 


2006- 


-04-18 


47 


3 


.2 


UGC 00238 


LINER 


00 


25 


03 


35 


+31 


20 


42. 


6 


6796 


2005- 


02-01 


35 


5 


6 


M31 


LINER 


00 


12 


11 


32 


+41 


16 


08. 


5 


-300 


2005- 


02-01 


13 


7 


.9 


IRAS 01189+2156 


Sy2 


01 


21 


44. 


27 


+22 


12 


35. 


1 


17478 


2005- 


-02-01 


35 


5 


.3 


UM319 


Sy2 


01 


23 


21 


18 


-01 


58 


36. 





4835 


2005- 


09-21 


56 


9 


.2 


UGC 00987 


Sy2 


01 


25 


31. 


16 


+32 


08 


11 


4 


4658 


2005- 


■01-28 


31 


5 


.2 


NGC 0660 


LINER 


01 


43 


01 


70 


+13 


38 


34. 





850 


2005- 


■02-02 


13 


8 


.7 


UGC 01282 


Sy2 


01 


49 


29. 


87 


+12 


30 


32. 


6 


5221 


2005- 


■02-02 


38 


6 


.5 


UGC 01479 


Sy2 


02 


00 


19 


06 


+24 


28 


25. 


3 


4927 


2005- 


■02-02 


36 


5 


7 


UGC 01757 


Sy2 


02 


17 


23 


05 


+38 


21 


19 


9 


5254 


2005- 


•01-28 


30 


5 


.5 


UGC 02456 


Sy2 


02 


59 


58. 


59 


+36 


49 


14. 


3 


3605 


2005- 


•02-01 


42 


5 


9 


UGC 02638 


LINER 


03 


17 


02 


21 


+01 


15 


17. 


9 


7123 


2005- 


09-21 


51 


10 


NGC 1320* 


Sy2 


03 


21 


18. 


70 


-03 


02 


32 


2 


2663 


2005- 


09-30 


17 


5 


.9 


MCG -02-09-040 


Sy2 


03 


25 


04. 


94 


-12 


18 


27. 


8 


4495 


2005- 


09-30 


48 


6 


.0 


SBS 0811+584 


Sy2 


08 


16 


01. 


30 


+58 


20 


01 





7584 


2005- 


09-30 


13 


5 


.6 


NGC 2683 


Sy2 


08 


52 


11 


12 


+33 


25 


13 


7 


411 


2005- 


09-30 


11 


6 


.5 


NGC 3190 


LINER 


10 


18 


05 


77 


+21 


49 


55. 


8 


1271 


2005- 


■10-01 


46 


5 


.9 


UGC 05613 


LINER 


10 


23 


32. 


50 


+52 


20 


30. 





9620 


2005- 


02-01 


34 


6 





NGC 3593 


Sy2 


11 


11 


37 


00 


+12 


19 


04 





628 


2005- 


■10-01 


17 


6 


.5 


NGC 3628 


LINER 


11 


20 


16. 


95 


+13 


35 


20. 


1 


843 


2005- 


02-01 


52 


8 


.1 


NGC 3753 


AGN 


11 


37 


53. 


90 


+21 


58 


53. 





8727 


2005- 


02-01 


40 


5 


.5 


MCG +03-38-017 


AGN 


11 


17 


53 


39 


+19 


01 


37 


1 


12506 


2005- 


02-01 


35 


6 


.0 


IRASF 15033+2617 


Sy2 


15 


05 


27. 


93 


+26 


05 


29 


3 


16579 


2005- 


02-01 


35 


5 


.9 


IRASF 15588+6522 


AGN 


15 


59 


18. 


95 


+65 


13 


57. 


8 


8975 


2005- 


■02-01 


31 


6 


.0 


IRAS 16288+3929* 


Sy2 


16 


30 


32 


66 


+39 


23 


03 


2 


9161 


2006- 


■04-10 


39 


2 


.8 


UGC 10593 


Sy2 


16 


52 


18. 


87 


+55 


51 


19 


8 


8739 


2005- 


09-21 


70 


14 


UGC 10695 


LINER 


17 


05 


05 


56 


+43 


02 


35. 


1 


8328 


2005- 


■09-28 


51 


6 


.6 


NGC 6503 


LINER 


17 


19 


27. 


11 


+70 


08 


39 


6 


60 


2005- 


■09-28 


56 


6 


.5 


MCG +11-22-046 


Sy2 


18 


22 


02 


20 


+66 


36 


37 


7 


4393 


2005- 


■09-28 


57 


7 


.0 


CGCG 341-006 


Sy2 


18 


45 


26. 


23 


+72 


11 


01 


7 


13880 


2005- 


•02-01 


35 


6 





NGC 7013 


LINER 


21 


03 


33 


31 


+29 


53 


19 


3 


779 


2005- 


■02-02 


38 


5 


.9 


IC 1417 


Sy2 


22 


00 


21. 


61 


-13 


08 


19 


1 


5345 


2005- 


09-30 


53 


6 


.1 


NGC 7331 


LINER 


22 


37 


04. 


09 


+34 


24 


56. 


3 


816 


2005- 


02-01 


42 


7 


.2 


UGC 12201 


LINER 


22 


19 


09. 


55 


+34 


59 


30. 


5 


5055 


2005- 


■02-02 


34 


1 


9 


NGC 7466 


Sy2 


23 


02 


03 


12 


+27 


03 


09 


5 


7493 


2005- 


02-01 


32 


5 


1 


UGC 12318 


Sy2 


23 


02 


08. 


01 


+25 


40 


13. 


7 


9533 


2005- 


02-01 


34 


5 


.5 


NGC 7549 


AGN 


23 


15 


17. 


23 


+19 


02 


30. 


1 


4736 


2005- 


■02-02 


35 


5 


.5 


UGC 12519 


Sy2 


23 


20 


02 


75 


+15 


57 


10. 


6 


4378 


2005- 


■02-01 


38 


5 


.6 



"Newly discovered H2O maser emission hosts are indicated by a ^. 

''Optical positions used for GBT pointing. Uncertainties are typically less than 2". Newly discovered maser emission 
in NGC 1320 and IRAS 16288+3929 lies within 1<t of the catalog positions for the two nuclei. 
NGC 1320: 02000 = 03 ,i 24 m 48?70±0?02;<52000 = -03°02'32"3 ± 0'.' 3 (maser; VLA); 

a 2000 = 03' l 24 m 48?70 ± 0?08; <5 2 ooo = -03°02'32;'2 ± l!'3 (nucleus; Iskrutskie et al1l2006l l ■ 
IRAS16288+3929: 02000 = 16 h 30 m 32?66 ± 0?02; <5 20 00 = 39°23'03?'l ± 0('3 (maser; VLA); 

a 2000 = 16 h 30 m 32f65 ± 0?02; 5 2 ooo = 39°23'03?'13 ± 0"5 (nucleus: lAbazaiian et aTll2005f) . 

c Systemic velocity (optical heliocentric) used in tuning the GBT observing bands. Velocities obtained from the 
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NED. When multiple measurements are available, those with the smallest listed uncertainties were chosen, in most 
cases. 

d Average system temperature. Temperatures (and RMS noise) during 2005 January and February integrations 
have been increased by a conservative 30% factor to correct for apparent miscalibration of the noise-source. The cor- 
rection factor was determined using a receiver temperature of 16.5 K, 6K contribution due to spillover and microwave 
background emission, and archival average forecast ground temperature and 22 GHz opacity. 

e RMS noise in a 24.4kHz spectral channel corrected for atmospheric opacity (i.e., referenced to outside the atmo- 
sphere), for the dependence of antenna gain on elevation, and for uncertainty in the noise source power level (see note 
d. Calibration uncertainties are discussed in Section 2. 
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Table 2. Deep Integrations on Known H2O Maser Sources 



Galaxy 


Type" 


Q?2000 i> 


^2000 b 


v c 
u sys 


Date 


rp d 

Jsys 


lcr e 






(hh mm ss) 


(dd mm ss) 


(km s- 1 ) 




(K) 


(mjy) 


NGC449 (Mrkl) 


Sy2 


01 16 07.25 


+33 05 22.4 


4780 


2006-04-19 


48 


2.4 


NGC 2979 


Sy2 


09 43 08.65 


-10 23 00.0 


2720 


2006-02-09 


37 


2.0 


NGC 3735 


Sy2 


11 35 57.30 


+70 32 08.1 


2696 


2006-02-09 


37 


3.9 


UGC4203 (Mrkl210) 


Sy2 


08 04 05.86 


+05 06 49.8 


4046 


2006-04-24 


11 


3.3 


Arp299 (NGC 3690) 


AGN 


11 28 32.20 


+58 33 44.0 


3121 


2006-04-27 


16 


2.9 


NGC 4293 


LINER 


12 21 12.82 


+18 22 57.4 


893 


2006-04-10 


38 


2.8 


NGC 0235A 


Sy2 


00 42 52.81 


-23 32 27.8 


6664 


2006-02-01 


58 


3.8 


NGC 2824 


AGN 


09 19 02.22 


+26 16 12.0 


2760 


2006-01-27 


35 


2.5 


IRASF 19370-0131 


Sy2 


19 39 38.91 


-01 24 33.2 


6188 


2006-04-19 


48 


3.7 



"Activity type from the NED, except for NGC0235A, which is from Kondratko et al. 2006a. 

ft GBT pointing positions adopted from Kondratko ct al. 2006a for NGC0235A, NGC 4293 (cr = Of'3), from 
Greenhill et al. 2003b for NGC 2824, NGC 2979, IRASF 19370-0131 (cr = 0!'2), and from the NED for NGC 449 
(o- = 0('28), UGC4203 (cr = Of 75), Arp299 (NGC 3690, cr = 10"), NGC 3735 (cr = K'3). 

c Systemic velocity (optical heliocentric) obtained from the NED used in tuning the GBT observing bands. 
Measurements for individ ual galaxies follow (in km s -1 ). NGC 449 (Mrk l): 4750 ± 16 (optical; Falco et al. 
1999), 4774 + 9 (optical: [Keell Il996h. 4823 ± 6_(HI; IBottinelli et al.lll990l) ; NGC2979: 2664 + 5 (HI; Saun- 

ioaD: 



ders et al. 2000), 2686 ± 8 (HI; IPaturel et alj [20031) , 2720 + 15 (optical; Fisher et al. 1995), 2731 ± 26 
(optical; iFairall et ai71ll992r); NGC 3735: 2597 ± 42 (optical; ISandagel Il978ri, 2684 + 39 (optical; Falco ct al. 



1999), 2696 + 6 (HI; Istavelev-Smith fc Daviesl|l988ri , 2696 ± 15 (HI; IFisher fc Tully|[l98ll) . 2738 ± 39 (optica l; 

" -9791 



Ide Vaucouleurs, de Vaucouleurs, fc Nietoll979n ; UGC 04 203 (Mrkl210): 3988 + 33 (optical; IFisher et al.lfl995h 
4046 ± 9 (HI; IBottinelli et al.lll990h , 4065 ± 42 (optical; Falco et al. 1999); Arp 29 9 (NGC 3690): see Hibbard fc 
Yun 1999 concerning nomenclature): 3040 + 27 (optical; Garcfa-Marfn et al.|[20o3 ). 3 080 + 40 (radio; IZhao et alj 
Il997l) ; Arp 299 (IC694): 31 20 ± 50 (radio; IZhao et alj|l997tl . 3121 + 27 (optical; [Garci'a-Marm et alj l2006h ; 



NGC 4293: 717 + 58 (optical; Ide Vaucouleurs et al.lll991), 867 ± 29 (opti cal; J. Huchra, priv ate communication) , 
883 + 54 (optical; Falco et al. 1999), 893 + 8 (optical; Idi Nella et al]|l995ri , 933 + 38 (optica l; iHuchra et al.| | l983 ); 
NGC0235A (ESQ 414-G016): 6692 + 37 (optical; da Costa et al. 1998), 6692 + 28 (optical; Ida Costa et al1l99lh 



6712 + 28 (optical; Wegner et al. 2003), 6765 ± 26 (optical; Ijones et aljliboih ; NGC 2824: 2760 + 23 (optical; 
Huchra et allll983h . 2755 ± 44 (optical; Falc o et al. 1999); IRASF 19370-0131: 6188 + 36 (optical; Istrauss et alj 



1992h . 6262 + 50 (optical; Ijones et al.1l2005ll 

d Average system temperature (not corrected for opacity). 



e RMS noise in a 24.4-kHz spectral channel corrected for atmospheric opacity (typically from 0.02 to 0.04) and 
for the dependence of antenna gain on elevation. 
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Table 3. Monitoring Log for Secular Drifts 



Source 


Epoch 


Time Monitor 11 


Mean Spectrum 6 


CG211 C 


2003- 


12- 


■08 


X 






2005- 


■01- 


■15 


X 


X 




2005- 


09- 


■28 


X 


X 




2006- 


■01- 


22 


X 


X 


NGC 6264 


2005- 


■01- 


■27 




X 




2005- 


10- 


■17 


X 


X 




2006- 


■02- 


■01 


X 


X 




2006- 


06- 


■05 


X 




VV340A 


2003- 


12- 


13 


X 






2005- 


■01- 


■15 


X 


X 




2005- 


06- 


■07 


X 






2005- 


10- 


■01 


X 


X 




2006- 


■02- 


■01 


X 


X 



"Epochs over which spectral features were tracked, enabling estimation of 
accelerations (Figures 6,7). 

^Epochs used to generate the weighted mean spectrum (Figure 5). For each 
source, the best three spectra were selected, with particular attention to the 
flatness of the spectral baseline. 



c The CG designation re fers to the Case low-dispersion spectroscopic survey 
jPesch fc Sanduleak|[l98i . and later papers by the same authors). 
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Table 4. Secular Drifts in Low- Velocity (Systemic) Emission 



Source 


Epoch" 




Velocity 


Drift 1 * 


Epochs 






(km s" 1 ) 


(km s _1 yr — 1 ) 




CG211 e 


2003 Dec 08 


1.7 


14940.92 ± 0.06 


0.405 ± 0.005 


1 








14943.2 ± 0.4 


0.90 ±0.10 










14951.9 ± 0.1 


0.93 ± 0.03 










14954.70 ± 0.06 


0.58 ±0.03 










14963.9 ± 0.1 


0.85 ±0.07 










14967.7 ±0.1 


1.55 ±0.04 










14969.9 ± 0.1 


1.88 ±0.07 










14976.0 ± 0.2 


0.11 ±0.05 










14949.2 ± 0.9 


1.45 ±0.09 






2003 Dec 08 


1.0 


15025.5 ±0.1 


0.65 ±0.04 












(avg. = 0.9)^ 




NGC6264 


2005 Oct 17 


1.2 


10164.0 ± 0.2 


1.9 ±0.2 


3 








10164.4 ± 0.3 


1.7 ±0.5 










10168.2 ± 0.3 


2.1 ±0.3 










10174.7 ±0.3 


1.7 ±0.6 










10181.9 ± 0.3 


2.1 ±0.6 










10184.7 ±0.4 


2.0 ±0.4 










10186.0 ± 0.1 


2.1 ±0.1 










10188.7 ±0.2 


2.0 ±0.2 










10191.0 ± 0.2 


1.2 ±0.4 












(avg. = 2)f 




VV340A e 


2003 Dec 13 


1.4 


10199.6 ± 0.1 


3.8 ±0.1 


5 



"Reference epoch for component velocities. 

''Reduced chi-squarcd for the least-squares decomposition of line profiles into Gaussian compo- 
nents, in which reference-epoch velocities, secular drifts, time- varying widths, and light curves are 
simultaneously fit for each Gaussian component using a time-series of spectra. Only components 
with signal-to-noise > 3 for two or more epochs are reported. 

c Reference epoch for fitted velocity. 

^Velocity drift of a Gaussian component. 

e Low- velocity plateaus were modeled by an FWHM ~ 25 km s" 1 stationary Gaussian compo- 
nent for CG211 and a first-order polynomial for VV340A. These components are not reported 
here. 

f Unweighted arithematic average for use in the calculation of entries in Table 5 (see discussion 
in Section 4.1). 
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Table 5. Estimate of 
Enclosed Mass and Disk Radii 



Galaxy 


Radius a 


Mass a 




(pc) 


(1O 7 M ) 


CG211 


0.4 


3 


NGC 6264 


0.2 


2 


VV340A 


0.3 


7 



a Because of the uncertainties due to 
the lack of VLBI data as discussed 
in Section 4.1, the accuracies of the 
masses and radii may be as large as half 
an order of magnitude. 
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Table 6. Inclinations of Galaxies that have Candidate or 
Confirmed Disk Masers 



Evidence for Disk Maser 6 
Source T b cos -1 (b/a) c i m od d Spectrum Map 

(°) (°) Ref. Ref. 



Nut 4945 1 


6 


70 


Circmus J 


3 


64 


NGU 4258 J 


1 


67 


NGC 1386^ 


-0.6 


68 


NGU 1068 J 


3 


32 


NGC 3079^ 


7 


80 


NGC 4388 


3 


77 


NGC 3735 


5 


78 


NGC 1320 


1 


70 


NGC 5728 


1 


51 


IC 2560-f 


3.3 


51 


NGC5793 


3 


70 


ngc 


1 


21 


NGC 591 





34 


Mrk 1419 


1 


46 


NGC 6926 


1 


16 


NGC 5495 


1.6 


37 


NGC 6323-f 


2 


71 


NGC 2979 


1 


51 


UGC 3789-'' 


2 


29 


UGC3193 


2.3 


68 


UGC 4203 


1 





NGC 449 


5 


52 


ESO 269-G012 


2 


39 


NGC 6264 




51 



86 


1 


1 


69 


2, 3 
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"Published disk maser candidates, in order of increasing recessional velocity. 
Missing isophotal sizes or merger activity exclude IRASF 22265-1826; CG211, 
VV340A, 3C403, NGC 17. 

''Morphological classification. 

inclination from axial ratio b/a. 



Modified inclination l lFouQue et al.lll990l) . Numbers in parentheses are ob- 
tained by adding lcr to T to obtain a real inclinations. 

e References for identifica tion of disk maser cand idates via spectroscopy : (1) 
Greenhill et al. (1997c); (2llNakai et a.1.1 dl995l"); OllGreenhill et al.l Jl997al') ; (41 
iNakai. Inoue. fc Mivoshil lll993ll; ( 5) iBraatz et al. I Jl996h; (6 ) IClaussen fc Lol 
Jl986h; (71 lHagiwara et alj i2002l1; (81 iBraatz et al.l d2004ll ; (9) this work 



(101 | Hagiwara et alj lll997 ); ( 111 IKondratko et al 
feOOal; (131 IGreenhill et alj ll2003rj) ; (141 iBraat: 

Kondratko et al. I l l2006bl l. References for confirmation of disk maser 



l 2006aT); ( 121 iHenkel etal] 
fc Gugliucci I 1 I2OO8II : (151 



iden- 



tification via VLBI study: 



IKond ratko et ah ~| 



(2005 



(161 IGreenhill et al. 1 J2003all; (171 


Argon et alj 


1 Jl997all: (191 Greenhill & Gwinnl 1 


1997)1; (201 


(211 llshihara et alj ||2001); (221 IKondratko et alj 



j2008ll ; (23l lBraatz et al~l i bOOTt) ; (24l lReid et alj J2009h 

^Disk maser status confirmed by VLBI observations. 



